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Abs  trac  t 


This  report  documents  the  results  of  the  experimental 
tests  performed  on  a  T-38  horizontal  stabilator  and  the 
investigation  of  a  two-dimensional  finite  element  model  for 
prediction  of  static  displacement  using  NASA  structural 
analysis  (NASTRAN)  computer  program.  The  finite  element  model 
was  used  in  previous  and  concurrent  theses  efforts  with  regard 
to  T-38  flutter  prediction  via  NASTRAN  modal  analysis  methods. 
It  will  also  be  used  in  future  analyses  by  a  Using  Air  Force 
Agency  for  flutter  prediction. 

The  experimental  test  set-up  for  the  stabilator  measured 
the  displacements  in  the  downward  (water  line)  direction  for 
ten  load  conditions  with  the  displacements  measured  at  twenty- 
five  locations  for  each  of  the  loading  conditions.  Precision 
load  cells  were  used  to  measure  the  applied  load  and  the 
resistance  to  rotation  at  the  actuator  arm.  The  boundary 
conditions  of  the  tests  represented  an  asymmetric  loading 
condition  for  the  T-38  assemblage. 

The  need  for  an  improved  two-dimensional  finite  element 
model  was  based  on  the  results  of  previous  flutter  investi¬ 
gations  using  the  same  two-dimensional  model.  The  model  was 
to  be  used  to  analyze  degradations  in  flutter  speed  due  to 
repair.  It  was  postulated  that  by  optimizing  the  model 
characteristics  for  the  static  displacement  cases  that  the 
improved  model  would  better  predict  the  flutter  speed  also. 


x 


The  results  of  the  experimental  tests  were  compared  with 
those  developed  using  NASTRAN.  Several  of  the  model  charac¬ 
teristics  were  altered  and  investigated  as  to  their  effect 
on  the  prediction  of  the  static  displacement.  The  analytical 
displacements  were  compared  with  those  achieved  experimentally. 

Finally,  development  of  an  optimization  procedure  is 
presented,  as  well  as  its  application  to  tuning  the  finite 
element  model.  The  development  of  the  optimization  routine 
was  based  on  current  information  for  optimization  techniques. 

A  comparison  is  drawn  between  the  error  value  used  herein 
and  current  optimization  theory. 
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INVESTIGATION  OF  A  TWO-DIMENSIONAL  MODEL 


FOR  THE  PREDICTION  OF  STATIC  DISPLACEMENT 
FOR  T-38  HORIZONTAL  STABILATORS  USING  NASTRAN 


I  Introduction 


Purpose 

The  horizontal  stabilator  of  the  United  States  Air  Force 
T-38  Talon  supersonic  jet  trainer  (henceforth,  referred  to 
as  stabilator)  is  the  subject  of  this  thesis.  San  Antonio 
Air  Logistics  Center  (SAALC)  currently  has  primary  engineer¬ 
ing  responsibility  for  repair  and  maintenance  of  the  T-38 
stabilator.  SAALC  has  the  requirement  for  an  improved  method 
of  examining  the  vibration  modes  and  frequencies  for  flutter 
speed  prediction  of  repaired  stabilators.  A  two-dimensional 
finite  element  model  (FEM)  was  developed  to  further  analyze 
the  modal  frequencies  of  the  stabilator  using  NASTRAN.  The 
use  of  modern  computer  numerical  analysis  methods  (NASTRAN) 
represented  an  improvement  over  traditional  methods  of  analy¬ 
sis  using  strip  theory.  However,  the  accuracy  of  NASTRAN' s 
capability  to  predict  static  displacements  and  modal  fre¬ 
quencies  (thus  predicting  the  flutter  speed)  is  directly 
related  to  the  input  information  contained  in  the  bulk  data 
deck  of  the  FEM.  Therefore,  it  is  the  purpose  of  this  thesis 
to  investigate  the  two-dimensional  finite  element  model  char¬ 
acteristics  for  predicting  static  displacements  for  the  T-38 
horizontal  stabilator  using  NASTRAN. 


Background 

Two  stabilators  were  built  by  Northrop  Aircraft,  Inc. 
(Northrop)  for  T-38  usage,  designated  Series  2  and  Series  3. 
Initial  development  and  testing  was  performed  on  the  Series  2 
stabilator  containing  three  auxiliary  ribs  which  were  deleted 
from  the  Series  3  design.  However,  the  stabilator  in  use  cur¬ 
rently  is  the  Series  3  design.  Thus,  a  void  in  the  data  for 
proper  deflection  analysis  based  on  the  Series  3  configuration 
existed  (for  both  static  and  dynamic  conditions).  The  research 
and  analysis  presented  herein  was  developed  using  the  Series  3 
stabilator  only  (Fig  1,  Table  I). 

The  two-dimensional  FEM  has  been  used  by  previous  thesis 
students  John  0.  Lassiter  (Ref  1)  and  Roger  K.  Thomson  (Ref  2) 
for  use  in  flutter  analysis.  In  each  of  the  above  studies, 
reference  was  made  to  the  discrepancy  in  static  displacements 
(from  NASTRAN  verification  computations)  when  compared  with 
published  Northrop  test  results.  Speculation  was  made  as  to 
the  probable  causes.  Cited  were  possible  FEM  errors  and  the 
lack  of  applicable  experimental  data  for  comparison  of  static 
deflection  verification.  It  was  a  logical  extension  of  the 
previous  theses  efforts  to  perform  the  needed  experimental 
tests  for  static  deflections  and  attempt  to  refine  the  FEM 
for  the  static  displacement  case. 

It  is  assumed  the  reader  has  an  adequate  background  or 
working  knowledge  of  NASTRAN  or  other  finite  element  methods 
for  structural  analysis  with  pax ^icular  emphasis  in  the 
solution  techniques  employed  in  static  deformation  problems. 
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Further,  it  would  be  of  benefit  in  evaluating  the  results  if 
the  reader  were  knowledgeable  in  the  area  of  flutter  analysis; 
inasmuch  as,  the  end  goal  of  the  FEM  is  to  predict  flutter 
speed  for  the  T-38  based  on  stabilator  vibration. 

Statement  of  the  Problem 


Investigate  the  geometric  and  physical  characteristics 
of  the  two-dimensional  FEM  used  with  NASTRAN  to  predict  static 
displacement  for  a  T-38  Series  3  stabilator,  with  the  goal  of 
minimizing  the  difference  between  analytical  and  experimental 
results . 

Approach  to  the  Problem 


As  inferred  earlier,  the  lack  of  experimental  data  for 
deflection  measurement  of  the  Series  3  stabilator  had  created 
a  problem  in  verifying  the  quality  of  the  two-dimensional 
finite  element  model.  The  initial  step  to  resolving  the 
problem  was  to  perform  the  necessary  experimental  tests  on 
the  Series  3  stabilator.  Once  the  tests  were  accomplished 
and  the  results  reduced,  the  next  step  was  to  accomplish  a 
comparison  between  the  results  produced  analytically  using 
the  Series  3  FEM  with  NASTRAN  and  the  results  achieved 
experimentally.  Lastly,  effects  of  varying  model  charac¬ 
teristics  (eg.  Young's  Modulus,  element  inertia,  thickness, 
etc.)  were  to  be  investigated  with  the  goal  of  developing 
an  optimum  FEM  for  the  NASTRAN  prediction  of  static  displace¬ 
ment.  The  optimized  FEM  would  then  be  used  for  flutter 
prediction . 
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II  Series  3  Static  Deflection  Tests 


Introduction 

The  experimental  tests  were  performed  on  a  Series  3 
stabilator  (Fig  1).  The  tests  paralleled  that  developed  by 
Northrop  (Ref  4)  for  the  Series  2  stabilator.  Significant 
deviations  from  the  Northrop  report  were  the  pad  size  (96  sq 
in  (Ref  4)  vs  49  sq  in  (used  herein)),  several  of  the  load 
conditions  (maximum  experimental  load  was  1200  lbs),  the  use 
of  a  load  cell  in  place  of  the  control  system  actuator  ram, 
and  the  boundary  conditions  (use  of  self-aligning  bearing  in 
two  places  in  a  rigid  structure).  Because  of  the  differences, 
it  was  impossible  to  correlate  the  two  tests  although  the  load 
locations  and  displacement  measurement  locations  were  the  same. 
The  applied  load  was  measured  in  twenty  percent  increments 
and  data  recorded  for  each  increment.  By  recording  the  data 
at  each  twenty  percent  of  load, it  was  possible  to  check  the 
linearity  of  any  particular  displacement  measured.  Also, 
each  load  was  cycled  a  minimum  of  two  times  to  reduce  any 
possible  hysteresis  of  the  system  before  recording  data. 
However,  only  the  maximum  load  condition  deflections  were 
averaged  and  used  for  analysis  comparison.  Therefore,  any 
further  reference  to  a  load  or  displacement  is  based  on  the 
maximum  load  value  and  its  corresponding  displacements. 
Presented  in  this  section  is  a  description  of  the  experi¬ 
mental  test  apparatuses,  data  acquisition  equipment, 
procedures  used  in  acquiring  the  data,  and  data  reduction. 
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Test  Apparatus 

The  experimental  tests  were  performed  using  the  resources 
of  the  Air  Force  Wright  Aeronautical  Laboratories  (AFWAL) , 
building  65,  Wright-Patterson  AFB ,  Ohio  (Table  II).  The  struc¬ 
ture  used  in  supporting  the  stabilator  for  the  tests  was  made 
of  steel  and  aluminum  components  (Fig  2).  "C"  channel  beams 

with  riveted  spacers  formed  "I"  beam  sections.  These  were 
bolted,  with  separate  corner  pieces,  to  form  the  support  struc¬ 
ture.  The  stabilator  was  supported  using  self-aligning  bearings 
in  aluminum  fixtures.  These  were  located  on  the  torque  tube 
at  HSS  25.07  and  HSS  0.00,  outboard  bearing  and  aircraft  center- 
line  positions,  respectively.  The  self-aligning  bearings  were 
the  same  as  those  used  in  the  aircraft.  This  represented  an 
asymmetric  loading  condition.  The  stabilator  was  mounted  such 
that  the  stabilator  planform  was  parallel  to  the  floor.  This 
allowed  the  displacements  to  be  measured  vertically. 

A  two  inch  diameter  ram  was  used  to  provide  the  displace¬ 
ment  loads.  The  Revere  precision  load  cell,  connected  to  the 
ram,  was  calibrated  to  1000  pounds  with  a  load  capacity  of 
1500  pounds.  The  pad  used  for  load  transfer  was  approximately 
fifty  square  inches  (7  in  by  7  in).  A  rod  end  joining  the 
load  cell  to  the  pad  allowed  the  pad  to  rotate  as  necessary  to 
meet  the  surface  angle  of  the  stabilator.  Thus,  deviation  of 
the  ram  from  the  plumbed  position  was  insignificant.  Hydrau¬ 
lic  pressurization  of  the  ram  provided  stabilator  displacement 
(downward).  However,  it  was  possible  to  begin  and  end  with 
zero  displacement  calibrations,  since  the  ram  could  be  lifted 
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Data  Acquisition 

System  Dayton  Scientific  RTP-7400  Sequential  send/receive 

Digital  Computer  Digital  Equipment  PDP-11  32K,  on-line  data 

_ _  Corporation  processing 


from  the  surface  of  the  stabilator  using  return  hydraulic 


pressure . 


Figure  2.  Experimental  Apparatus  View  1 

Measurement  of  the  displacement  was  made  using  potentio- 
metric  displacement  transducers.  The  range  of  displacement 
made  it  necessary  to  use  three  sizes  of  transducers  (1,  2, 
and  3  in).  Twenty-five  transducers  were  attached  to  a  sub- 
frame  located  beneath  the  stabilator  (Fig  3).  Each  transducer 
position  was  found  using  a  plumb  line  from  the  attachment 
point  on  the  stabilator.  The  transducers  were  attached  to 
the  stabilator  via  .016  inch  stainless  steel  transmission  wire. 
The  twenty-five  positions  corresponded  to  those  in  Reference  4 
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testing.  A  manual  jack  located  beneath  the  outboard  bearing 
support  structure  served  to  enforce  the  zero  bearing  dis¬ 
placement  boundary  condition.  Force  in  the  downward  direction 
(to  enforce  zero  vertical  displacement)  was  applied  to  the 
centerline  bearing  by  way  of  tightening  adjustment  bolts  on 
the  support  structure. 


**&.**»*  -  *-f  -A.'  . 


Figure  3.  Experimental  Test  Apparatus  View  2 

Movement  of  the  bearings  (outboard  and  centerline)  was 
monitored  for  vertical  displacement  via  two  inch  dial  indica 
tor  gages  which  read  to  one  thousandth  of  an  inch.  The 
support  brackets  for  the  gages  were  attached  to  an  assumed 
rigid  structure  (Fig  4).  A  Toroid  precision  load  cell  was 


used  to  measure  the  reaction  loads  at  the  control  actuator 
arm.  The  spring  constant  for  the  load  cell  was  5.0  x  10^ 
pounds  per  inch.  It  was  calibrated  to  2500  pounds  with  a 
load  capacity  of  3750  pounds. 


Figure  4.  Torque  Tube  Boundary  Condition 
Measurement  Apparatus 

Hydraulic  pressure  was  obtained  using  manual  pumps  with 
internal  reservoirs  (Fig  5).  The  assembly  was  manufactured 
by  AFWAL.  Pressure  gages  monitored  the  line  pressure  as  a 
safety  precaution.  Fine  adjustment  of  the  hydraulic  pressure 
was  obtained  using  a  screw  device  to  create  small  changes  in 
hydraulic  line  pressure. 
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Figure  5.  Hydraulic  Pressure  Source 
and  Monitor  Screen 

Data  Acquisition  Equipment. 

A  video  screen  was  used  to  monitor  the  load  on  the 
stabilator  at  the  hydraulic  pump  station.  As  a  result  of 
on-line  data  processing  (in  real  time)  it  was  possible  to 
monitor  the  load  (measured  by  the  load  cell),  several  dis¬ 
placements  and  adjust  line  pressure  to  the  ram  at  the  same 
station  (Fig  5). 

A  separate  grouping  of  apparatuses,  shown  in  Figure  6, 
was  assembled  for  data  acquisition.  The  TEC  terminal / screen 
was  used  to  display  the  load  cell  values,  displacements,  and 
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interface  with  the  PDP-11  computer  for  data  processing.  The 
PDP-11  digital  computer  was  located  in  the  same  test  building. 
A  sample  rate  of  five  samples  per  second  was  selected  for  all 
tests,  in  as  much  as  the  tests  were  static  displacement  tests. 
Separate  power  and  balance  modules  were  used  for  each  load 
cell  and  each  transducer.  An  initial  no  load  displacement 
reading  of  each  transducer  prior  to  load  application  was  used 
as  a  zero  reference.  Zero  displacement  setting  for  each  trans¬ 
ducer  was  accomplished  by  displaying  the  voltages  sent  to  the 
PDP-11  digital  computer  on  the  real  time  peripheral  unit  and 
adjusting  the  appropriate  module  for  zero  displacement  (out¬ 
put  voltage).  A  separate  digital  voltmeter  was  used  for 
calibration  readings  for  the  load  cells. 


Figure  6 


Data  Acquisition  Console 
13 
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The  Digital  Equipment  Corporation  PDP-11  digital  com¬ 
puter,  housed  in  building  65,  was  used  for  data  storage, 
sampling,  and  processing.  The  unit  allowed  for  data  to  be 
read  on  the  video  monitors  prior  to  its  transfer  to  tape 
storage  and  post-processing.  The  data  was  processed  on-line 
in  engineering  units  in  real  time. 

Experimental  Procedure 

The  stabilator  was  placed  in  the  support  structure,  pre¬ 
viously  discussed,  using  self-aligning  bearings.  The  mid-thick¬ 
ness  plane  (planform)  of  the  stabilator  was  parallel  to  the 
floor.  The  upper  and  lower  surface  load  application  and  dis¬ 
placement  locations  were  based  on  dimensions  contained  in 
Reference  3. 

The  position  of  each  transducer  was  determined  using  a 
plumb  line  from  the  underside  attachment  fixtures.  The  dis¬ 
placement  transducers  were  connected  to  the  fixtures  on  the 
stabilator  by  fine  stainless  steel  wire.  The  position  of  the 
torque  tube  height  was  checked  and  the  dial  indicators  zeroed. 
The  twenty-five  transducers  were  then  zeroed  (relative  position 
of  the  variable  resistance  transducer  as  measured  by  the  real 
time  peripheral  unit  was  recorded  by  the  PDP-11  as  the  start 
point  for  displacements  with  no  load).  The  bridge  balance  on 
the  load  cells  were  checked  and  calibrated  similarly. 

The  actuator  position  was  determined  by  a  plumb  line  from 
the  overhead  support  beam  to  the  position  previously  marked 
for  a  particular  load  application.  Cross  lines  marking  the 
position  were  used  to  center  the  pad.  Gradual  hydraulic 
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pressure  was  applied  until  contact  of  the  rubber  pad  was  made 
on  the  stabilator.  Hydraulic  pressure  was  increased  and  the 
applied  load  monitored  on  the  video  screens.  The  displace¬ 
ment  of  the  torque  tube  was  monitored  and  forced  to  zero 
vertically,  using  the  mechanical  jack  and  centerline  bearing 
support  adjustment  fixture.  The  data  was  recorded  for  the 
applied  load,  reaction  load  at  the  actuator  arm,  and  the 
twenty-five  displacements  for  the  particular  load  case.  The 
load  was  then  removed  gradually,  with  the  torque  tube  dis¬ 
placement  continuously  being  held  at  zero  vertically.  The 
zero  load  residual  displacements  for  the  transducers  were 
displayed  (hysteresis  effects).  The  process  of  calibrate, 
load,  measure,  unload  was  then  repeated.  When  the  system 
reached  a  hysteresis  effect  of  less  than  five  thousandths  of 
an  inch  (for  the  worst  case),  the  data  was  saved  for  post¬ 
processing,  using  commands  from  the  TEC  unit.  One  additional 
run  was  made  for  repeatability  demonstration.  This  procedure 
was  repeated  for  ten  distinct  load  positions. 

Data  Reduction 

Data  recorded  from  the  experimental  tests  were  post- 
processed  into  tabular  listings.  The  recorded  data  was 
saved  at  a  rate  of  five  samples  per  second.  Since  the  tests 
were  static,  any  variation  in  displacement  was  accounted  for 
by  averaging  each  transducer's  output  to  establish  a  single 
value  for  displacement  for  a  particular  load  condition, 

A  separate  FORTRAN  program  was  used  to  expand  the  exper¬ 
imental  data  by  linear  interpolation  into  a  grid  arrangement 
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which  was  the  same  as  that  used  for  NASTRAN  computations. 

The  interpolated  displacement  results  were  then  used  as 
input  to  the  graphics  computer  code  GCSNAST. 

Results 

The  experimental  tests  were  successful  with  respect  to 
providing  data  for  static  deflections  of  the  Series  3  stabi- 
lator.  A  correction  to  the  input  load  measured  at  the  ram 
was  necessary,  since  the  bridge  balance  for  the  load  cell 
was  accomplished  with  the  pad  suspended  on  the  end.  With 
this  correction  applied,  good  correlation  was  found  between 
the  applied  load  and  the  measured  reaction  load  at  the 
actuator  arm. 

Results  of  initial  trial  tests  indicated  a  need  for  fur¬ 
ther  support  on  the  torque  tube  bearings  to  enforce  zero 
vertical  displacements.  The  use  of  a  jack  under  the  outboard 
bearing  and  a  support  structure  with  adjusting  bolts  above 
the  centerline  bearing  solved  the  problem. 

The  measured  displacements  were  averaged  for  each  load 
condition.  The  results  are  presented  in  Appendix  A.  The 
displacements  were  then  linearly  interpolated  to  an  expanded 
set  of  displacements.  The  linearly  expanded  set  was  then  used 
with  the  computer  graphics  code  GCSNAST  to  produce  figures 
(App  B).  The  figures  illustrate  the  deformed  shape  and  were 
useful  in  comparing  the  analytical  results  (which  are  dis¬ 
cussed  in  the  following  chapter). 
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Ill  Stabilator  Finite  Element  Model  Analysis 


Introduction 

The  two-dimensional  finite  element  model  had  been  used 
in  previous  theses  by  Lassiter  (Ref  1)  and  Thomson  (Ref  2) 
as  mentioned  in  Section  I.  A  detailed  description  of  the 
initial  development  of  the  FEM  can  be  found  in  Reference  1, 
Section  II  and  Appendix  A.  Those  items  significant  to  the 
development  of  this  analysis  will  be  presented  here  also. 

The  data  necessary  to  describe  the  physical  characteristics 
of  the  stabilator  for  NASTRAN  (part  of  the  Bulk  Data  Deck) 
was  generated  from  a  separate  FORTRAN  program  (BDATA).  A 
description  of  the  data  generating  program  is  presented  in 
Appendix  C.  The  elemental  substructures  used  to  model  the 
stabilator  were  triangular  membrane  and  bending  elements 
(CTRIA1),  quadrilateral  membrane  and  bending  elements  (CQUAD1), 
and  simple  beam  elements  (CBAR)  (Ref  5).  In  the  interest  of 
clarity  and  congruity,  the  FEM  used  herein  was  referred  to  as 
the  influence  coefficient  model  in  the  above  theses.  The 
number  of  elements,  nodes,  and  physical  dimensions  for  base 
line  data  were  the  same  as  those  used  in  the  influence  coef¬ 
ficient  model  of  previous  theses. 

The  FEM  boundary  conditions  were  based  on  those  modelling 
the  experimental  tests  (not  the  same  as  the  actual  aircraft, 
previous  Northrop  tests,  or  other  theses  efforts). 
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Development  of  the  Finite  Element  Model 


Stabilator  Description.  The  T-38  Talon  stabilator  is 
a  lightweight,  honeycomb  structure  with  aluminum  panels 
(upper  and  lower  surface)  bonded  to  the  honeycomb  core  by 
epoxy.  The  aluminum  panels  are  riveted  to  substructures, 
such  as,  the  root  rib,  main  spar,  auxiliary  spar,  tip  rib, 
and  leading  edge  extrusion.  When  the  stabilator  is  installed, 
the  planform  tilts  four  degrees  down  from  the  root  rib  to  the 
tip  rib. 

Substructure  Elements.  NASTRAN,  utilized  in  solving 
aerodynamic  problems,  had  available  the  necessary  elements 
to  adequately  model  the  stabilator.  The  plate  bending 
elements  with  in-plane  stiffness  were  used  to  model  the 
skin-honeycomb  core  combination.  These  were  the  CQUAD1  and 
CTRIA1  elements.  Bar  elements  (CBAR)  were  used  to  model  the 
hinge  fitting  (torque  tube  and  main  spar),  the  auxiliary 
spar,  tip  rib,  leading  edge  forging,  root  rib  sections,  and 
trailing  edge  closure  (refer  to  Fig  1).  It  should  be  pointed 
out  this  combination  of  elements  forms  a  two-dimensional  FEM. 

The  influence  coefficient  model  was  developed  based  on 
15  quadrilateral  elements  spanwise  and  8  quadrilateral  ele¬ 
ments  chordwise.  Several  quadrilateral  elements  were  divided 
into  triangular  elements  (Fig  7).  The  physical  properties  were 
found  for  the  quadrilateral  elements  prior  to  being  divided  into 
triangular  elements.  The  values  of  the  physical  properties  at 
the  four  nodes  were  developed  using  the  elemental  properties  at 


the  center  of  the  element  (found  using  a  separate  algorithm). 
The  quadrilateral  elements  were  then  divided  into  triangular 
elements  as  required  to  form  the  influence  coefficient  mesh. 
Some  nodes  were  moved  to  prevent  poor  aspect  ratio  triangles 
or  for  location  of  load  application  purposes.  The  resulting 
triangular  elements  or  distorted  quadrilateral  elements  used 
the  same  physical  and  material  properties  as  were  developed 
at  the  center  of  the  quadrilateral  element  replaced. 

Elemental  properties  for  the  bar  elements  were  similarly 
formulated.  The  physical  and  material  properties  were  de¬ 
termined  at  the  center  of  the  element  and  used  for  the  entire 
element.  The  torque  tube  (a  part  of  the  hinge  fitting  -  Fig  1) 
was  developed  separately.  The  location  of  the  grid  points 
were  based  on  constraint  requirements  (eg.  bearing  locations, 
actuator  arm  location)  and  not  on  cross-sectional  properties. 
The  values  used  for  physical  properties  were  based  on  the 
dimensional  characteristics  of  the  section  (three  sections 
were  utilized).  Because  of  the  complexity  of  the  changing 
geometry  of  the  torque  tube  between  the  root  rib  and  the 
actuator  arm  location,  averaged  values  were  used  for  the 
sectional  properties. 

The  program  BDATA  generated  much  of  the  data  describing 
the  physical  and  material  properties  used  in  describing  the 
substructure  elements.  These  included: 

1.  Property  Cards  (PBAR,  PQUAD1,  PTRIA1)  identified  the 
physical  characteristics  of  the  element,  such  as  thickness, 
material  identifiers,  inertias,  mass. 
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2.  Material  Cards  (MAT1,  MAT2)  described  the  material 


characteristics,  such  as  Young's  modulus,  shear  modulus, 
Poisson's  ratio.  MAT1  was  used  for  isotropic  material. 

MAT2  was  used  for  orthotropic  material. 

3.  Grid  points  (GRID)  were  assigned  unique  geometric 
locations  in  terms  of  fuselage  station  (FS),  horizontal 
stabilator  station  (HSS)  and  water  line  (WL)  locations. 

GRID  also  identified  the  appropriate  coordinate  system  and 
unique  constraints  for  the  grid  point.  The  influence  coef¬ 
ficient  FEM  used  166  grid  points  (App  E) . 

4.  Due  to  the  four  degree  tilt  of  the  planform  in  the 
aircraft  installed  configurations,  a  second  coordinate  sys¬ 
tem  (C0RD2R)  was  used  for  the  planform  (Fig  8).  This  per¬ 
mitted  the  loads  to  be  applied  normal  to  the  planform  and 
the  displacements  analyzed  normal  to  the  undeformed  planform 
duplicating  the  loading  and  measured  displacement  angles  of 
the  experimental  tests. 

Significant  items  of  the  FEM  necessary  for  NASTRAN 
analysis  which  were  not  generated  using  BDATA  are  listed 
below. 

1.  The  Executive  Control  Deck  identified  the  solution 
format,  as  well  as  alterations  for  calculating  an  error 
value  (App  F) . 

2.  The  Case  Control  Deck  contained  information  relative 
to  controlling  subcase  solution  sets  and  output  requirements 

3.  The  single  point  constraints  identifying  constrained 
degrees  of  freedom  at  a  grid  point  were  used  for  fixed  boun¬ 
dary  conditions. 
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Figure  8.  Finite  Element  Model  Coordinate  Systems 


4.  A  scalar  elastic  element  was  used  to  model  the  load 
cell  actuator  arm. 

5.  The  experimental  values  of  the  twenty-five  node  dis¬ 
placements  were  input  for  use  in  the  Executive  Control  Deck 
Alter  routine  (App  F). 

6.  Several  other  matrices  were  input  for  data  extra¬ 
polation  and  matrix  partitioning  to  be  used  by  the  Executive 
Control  Alter  routine  (App  F) . 

7.  Lastly,  the  forces  were  input  with  directional 
cosines  and  appropriate  coordinate  system  identification. 

Boundary  Conditions.  The  boundary  conditions  imposed  on 
the  FEM  were  based  on  the  experimental  tests.  These  repre¬ 
sented  an  asymmetric  loading  condition  as  mentioned  in  the 
Introduction.  A  spring  constant  for  the  load  cell-actuator 
arm  was  calculated  from  available  data.  The  boundary  con¬ 
ditions  on  the  torque  tube  were  rigidly  fixed  using  single 
point  constraints  (SPC)  for  the  x  and  z  displacements  at  the 
outboard  bearing  and  x,  y  and  z  displacements  at  the  center- 
line  bearing  location. 

Static  Analysis 

For  static  analysis  of  the  FEM  and  the  ten  load  subcases 
(Table  III),  Rigid  Format  1  of  NASTRAN  was  used.  The  mathe¬ 
matical  formulation  of  the  solution  process  is  based  on  the 
linear  theory  of  elasticity.  The  FEM  stiffness  matrix  was 
formulated  via  the  Bulk  Data  Deck  (BDATA  generated  data  and 
hand  input  data  combined).  The  stiffness  matrix  was  then 
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TABLE  III 


¥ 


partitioned  with  respect  to  constrained  degrees  of  freedom 
and  free  degrees  of  freedom  to  be  solved.  Those  equations 
not  constrained,  form  a  subset  identified  as  the  1-set, 
which  is  the  solution  set  for  static  analysis  (Ref  6). 

The  influence  coefficient  model  was  developed  as  a  result 
of  the  requirement  that  the  load  (force)  be  concentrated  at 
a  grid  point  for  static  analysis.  The  set  of  forces  is  re¬ 
duced  to  a  1-set  for  nodal  forces  from  constraint  partitioning 
corresponding  to  the  stiffness  reduction.  This  reduces  the 
solution  process  to  that  of  satisfying  equilibrium  equations 
for  stiffness  versus  force.  The  resulting  equation  can  be 
expressed  as 

{P}  ='  [s]  {5}. 

The  vectors  P  and  5  represent  the  nodal  forces  and  displace¬ 
ments,  respectively,  and  the  matrix  S  represents  the  FEM 
stiffness  matrix.  For  solution  of  the  equations,  NASTRAN 
reduced  the  stiffness  matrix  into  its  upper  and  lower  tri¬ 
angular  factors.  Then  forward-backward  substitution  is 
performed  for  all  load  cases  having  the  same  constraints, 
which  was  the  case  here.  The  displacement  solution  for 
each  load  case  contained  three  translational  and  three 
rotational  values  for  each  grid  point  (except  those  dis¬ 
placements  constrained).  The  vertical  displacements  for 
twenty-five  grid  points,  to  be  compared  with  the  experimental 
displacements,  were  extracted  using  an  Alter  routine  in  the 
Executive  Control  Deck  section  of  the  Bulk  Data  Deck  (App  F) . 
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The  results  of  the  static  analysis  displacements  and  the 
experimental  displacements  were  compared. 

A  comparison  of  the  experimental  reaction  load  and 
computed  reaction  load  at  the  actuator  arm  was  made.  The 
reaction  loads  were  found  to  be  in  excellent  agreement 
after  the  magnitude  of  the  applied  load  was  corrected  by 
subtracting  the  weight  of  the  pad  suspended  from  the  load 
cell. 

Application  of  Optimization  Theory.  As  an  introduction 
to  the  development  and  implementation  of  an  error  function 
for  optimizing  the  difference  between  the  computed  displace¬ 
ments  and  the  experimental  displacements,  several  key 
optimization  principles  are  reviewed  (Ref  7).  Several 
current  computer  programs  used  for  implementation  of 
optimization  functions  are  also  reviewed  for  possible 
future  utilization. 

Similar  to  most  optimization  formulations,  the  error 
in  the  displacements  (or  effectiveness  of  the  system)  was 
calculated  from  a  number  of  individual  subcases  (sometimes 
only  one  case  is  optimized).  Each  of  the  subcases  may 
also  have  effectiveness  measurements  associated  with  them. 
For  the  problem  herein,  the  overall  effectiveness  might  be 
expressed  as 

n 

SE  =  2  z  -  (a,x) 

i=l  1  ~  “ 

where  z^  represents  a  performance  indicator  for  the  subcases 


and  is  a  function  of  a  and  x,  the  decision  and  state  vari¬ 
ables,  respectively.  This  relationship  is  used  to  find  the 


error  (SE  -  system  effectiveness).  In  this  case,  the  mini¬ 
mum  value  is  sought.  Each  of  the  subcases  can  be  expressed 
functionally  by  relationships  of  the  decision  variables, 
such  as 


Thus  , 


zi  =  f i Ca , x) 

zL  =  fj.CE,  I,  G,  J,  t,  etc) 


The  functions  developed  for  subcase  performance  identifiers 
relate  the  effects  or  sensitivity  of  the  decision  variables 
(a)  and  state  variables  (x)  to  the  effectiveness  measure 
for  that  particular  subcase  (z^)  -  an  extensive  effort  for 
ten  subcases.  The  combination  of  decision  variables  and 
state  variables  must  be  capable  of  describing  the  state  of 
the  system  (the  error  for  the  particular  subcase  z^)  for  any 
variable  (a  and/or  x)  change,  which  in  turn  is  used  to  es¬ 
tablish  the  value  of  the  system  effectiveness. 

For  the  analysis  to  be  manageable,  it  is  desirable  that 

the  performance  index  vector  Z  (where  Z  =  { z ^ , • • • , zn) ) 

have  several  basic  properties.  The  vector  Z  must  be  complete 

< 

with  respect  to  the  system  effectiveness.  'This  reiterates 
the  necessity  for  the  decision  variables  and  the  state  vari¬ 
ables  to  describe  the  state  of  the  system  for  any  variable 
change.  should  convey  meaningful  results  (as  it  did 
herein  -  summing  all  the  inches  differences  in  displacements). 
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Z  should  be  decomposable  into  generally  independent  elements. 
This  relates  to  one  of  the  problems  found  in  the  analysis; 
that  is,  changing  one  decision  variable  and  having  z^  values 
change  with  conflicting  results  (some  improved  -  others  be¬ 
came  worse).  To  avoid  having  the  SE  biased,  the  relationships 
of  the  subcase  identifiers  (z^)  must  be  non-redundant .  This 
avoids  being  rewarded  (or  penalized)  twice  or  more  for  the 
same  variable  change.  The  size  of  the  n-vector  Z  should  be 
kept  to  a  minimum.  The  recommendation  here  is  to  pick  one 
dominant  in  spanwise  bending,  subcase  (2);  two  dominant  in 
chordwise  bending,  subcases  (8  and  10,  one  for  each  rotation 
direction);  one  mixed  mode,  subcase  (4);  and  one  close  to 
the  boundary,  subcase  (9).  This  requires  deleting  five 
subcases  (1,  3,  5,  6,  and  7).  This  should  make  the  optimi¬ 
zation  tractable,  as  well  as  feasible. 

Once  the  necessary  state  variables  are  established,  which 
accurately  define  the  state  of  each  subcase  performance  indica¬ 
tor,  the  value  of  the  decision  variables  can  be  bounded  for 
the  decision  space.  The  constraining  equations  (C  )  further 

O 

limit  the  number  of  iterations  necessary  for  a  numerical 
solution.  However,  they  should  be  of  sufficient  size  to 
adequately  consider  all  combinations  for  arriving  at  the 
optimum  SE  solution  set.  As  an  example,  considering  the  dev¬ 
elopment  of  an  equivalent  thin  rectangular  bar  in  rotation 
for  replacing  a  torsion  cell,  Young's  modulus  might  be  made 
to  vary  from  E  to  4E .  Expressing  this  as  an  inequality,  the 
relationship  would  be 


There  are  no  limits  placed  on  the  number  of  constraint 

equations  (C  )  which  can  be  utilized.  For  the  problem 
6 

herein,  it  is  recommended  the  range  be  as  constraining 
as  is  feasible,  or  extended  in  small  increments  as 
necessary . 

In  the  event  it  is  desirable  to  evaluate  the  problem 
on  the  basis  of  a  preferred  variable,  then  a  resulting 
hierarchy  (weighting)  of  the  variables  leads  to  a  complete 
equivalent  ordering  of  the  performance  indicators  (z^). 
Mathematically,  this  can  be  expressed  as  an  alternative 
(individual  realizations  of  the  performance  indicators) 
and  is  represented  by 

z1  y  Z2 

to  be  read  as  alternative  Z^1  is  preferred  to  alternative 
Z2  where  Z1  represents  the  set  of  performance  indices 
associated  with  the  ith  alternative.  Further,  this  implies 
the  functional  relationship  f(z^)  is  greater  than  or  equal 

2_i  2 

to  f(z^  )  for  all  l,  with  f(z.  )  strictly  greater  than  f(z^  ) 
for  at  least  one  i. 

Definitive  relationships  between  decision  variables 
(often  identified  as  state  equations)  further  decrease  the 
number  of  actual  independent  variables  that  can  be  used  for 
optimization.  These  are  represented  as  equalities  (ST^). 

One  such  example  is 


G 


0 


STh(E) 


E 

2(l+v) 


Only  one  variable  (E  or  G)  is  necessary,  the  other  is  redun¬ 
dant  (assuming  isotropic  material  properties).  The  use  of 
orthotropic  material  properties  would  allow  the  use  of  in¬ 
dependent  bending  (two  directions)  and  shear  moduli,  thus 
having  three  decision  variables  to  use  for  optimization 
(assuming  the  two  bending  moduli  are  not  coupled  by  Poisson's 
ratio)  . 

The  method  of  solution  outlined  from  Reference  7  is  opti¬ 
mization  using  multiple  objective  optimization  theory  (MOOT). 
Mathematically  the  MOOT  process  can  be  represented  as: 


Extremize  SE  =  extremize  f(Z(x))  =  f(extremize  Z(x)) 

subject  to:  Cg(a,x)  <0  g  =  1,2,... ,k 

STh(a,x)  =0  h  =  l,2,...,m 

Where  extremize  means  to  minimize  the  error  value  for  the 
problem  herein.  The  result  of  the  optimization  is  a  set  of 
variables  (a,x's)  which  form  a  non-dominated  solution  (%*) 
set  (NDSS)  (provided  a  solution  exists).  A  non-dominated 
solution  (NDS)  is  defined  as: 

is  a  NDS  iff  ^  (a,x)  C  X  3 

zr(£>2£)  >  zr(£*)  f°r  some  r  =  l,2,...,n 

zs(a,x)  £  ZS(Z*)  v  s  f  r 

Therefore,  y*  is  a  NDS  containing  optimum  values  for  both 
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the  range  space  of  decision  and  state  variables  (X);  such 
that  zr(a,x)  is  better  than  for  some  r  =  l,2,...,n 

decision  variable;  and  zs(a,x)  better  than  or  equal  to 
z  <£*)  for  all  s  not  equal  to  r.  Allowing  X  to  be  the  range 
space  for  all  possible  realizations  of  a  and  x  (from  possible 
combinations  of  r  values),  then  is  a  specific  realization 
from  within  the  range  space  and  is  a  candidate  solution  set. 

The  software  for  solving  the  equations  to  arrive  at  a 
solution  for  the  NDSS  exists  in  several  formats  depending  on 
the  type  of  optimization  to  be  performed.  Examples  of  such 
programs  are  PROCES  (Ref  8)  and  AESOP  (Ref  7)  developed  at  the 
Air  Force  Institute  of  Technology.  The  standard  format  of  the 
input  for  these  programs  are  of  the  following  form 

Extremize  Z(a,x)  =  WW  *  [»]{*} 

subject  to  C  (a,x)  <  0  for  g  =  l,2,...,k 

O 

and  ST^(a,x)  =  0  for  h  =  1,2 . . 

where  Z(a,x)  is  the  performance  vector,  a  and  x  are  the 

decision  and  state  variables,  respectively,  and  A  and  B 

are  matrices  of  coefficients  relating  the  decision  and 

state  variables  to  the  performance  indices.  The  solutions 

satisfying  the  above  relations  are  then  compared  with  one 

another  through  a  separate  algorithm  which  checks  for 

2  . 

dominance.  If  the  second  candidate  set  (a,x)  is  better 
than  the  first  candidate  set  (a,x)^,  the  second  is  retained 


t! 


as  a  new  NDS  (^*) .  This  process  continues  until  all  can¬ 
didate  solution  sets  have  been  compared.  The  result  is  a 
NDSS,  or  the  best  solution  set  based  on  the  relationships 
input,  constraint  equations,  and  state  equations.  Implemen¬ 
tation  of  the  vector  optimization  process  is  accomplished 
by  the  mathematical  formulation 

Extremize  SE  =  extremize  f(Z(a,x))  =  f(extremize  Z(a,x)) 


reposed  as 


Extremize  zi(aJx) 


Cg(£,x)  <  0 

ST,  (a.,x)  =  0 


subject  to: 


for  all  g 


for  all  h 


Formulating  the  solution  using  the  constraint  method, 
requires  all  other  performance  indices  to  be  expressed 
as  state  equations,  such  as 


zk(a,x)  "  rk  =  0 


for  all  k  f  l  above 


The  constraint  solution  process  optimizes  one  functional 
relationship  while  all  other  relationships  are  held  constant. 
Where  the  range  of  all  possible  realizations  for  rk  depend 
upon  the  range  of  the  individual  objective  functions. 
Optimization  is  accomplished  sequentially  for  all  possible 
values  of  r^  (and  combinations  thereof).  The  candidate 
solutions  are  compared  as  before  to  determine  if  a  NDS 
exists  (or  no  solution  set  is  found  due  to  the  problem 
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being  over  constrained).  Then,  the  system  effectiveness 
is  calculated  based  on  the  summation  of  the  z^(a,x)  func¬ 
tionals  from  the  optimum  solution  set  of  variables  (a,x). 

If  the  constraint  method  proves  to  be  ineffective  for  the 
problem  at  hand,  the  weighting  method  may  be  a  feasible 
alternative.  However,  the  weighting  method  may  arrive  at 
a  solution  set  at  the  cost  of  some  numerical  accuracy. 

Several  factors  deemed  the  above  programs  unusable  for 
application  (at  the  time  analysis)  to  the  problem  pre¬ 
sented  herein.  These  included: 

1.  Assuming  non-linear  relationships  of  the  decision 
variables  could  be  derived,  it  would  have  required  an  ex¬ 
cessive  number  of  relationships  to  develop  the  minimumization 
optimizing  functions  for  the  problem  herein.  None  of  the 
relationships  could  be  deduced  based  on  physical  properties. 

The  relationships  developed  would  have  to  relate  a  change  of 
decision  or  state  variable  with  a  displacement  difference 
(versus  pure  displacements  found  in  a  static  analysis  problem). 

2.  The  constraining  equations  necessary  to  search  the 
entire  domain  of  possible  realizations  were  unknown.  Thus, 
the  range  would  be  extensive  for  searching  for  the  optimum 
solution  set. 

3.  The  state  variable  relationships  were  also  unknown 
and  could  not  be  determined  from  physical  insi6ht  of  the 
problem. 

4.  AESOP  (ability  to  handle  non-linear  relationships) 
was  being  developed  concurrent  with  this  thesis  effort. 


However,  much  of  the  optimization  material  presented  was 
applied  in  deriving  a  method  of  error  analysis  for  the 
stabilator.  The  major  differences  are  noted  as  development 
of  the  error  value  formulation  is  presented. 


Error  Value  Formulation.  An  error  value 
on  optimization  principles  was  applied  to  the 
The  error  value  was  formulated  mathematically 
minimization  of  the  system  effectiveness  value 
the  minimum  value  of  each  of  the  subcases  (zj) 
the  overall  error  value  (E)  being  minimized. 


function  based 
problem  herein 
based  on  the 
In  theory, 
resulted  in 


E 


iq 

10 

minimum 

z  . 

= 

2  minimum  z . 

J 

j=l  J 

Where  E  is  the  error  (equivalent  to  SE  in  this  application) 
in  the  FEM  (analytical  versus  experimental).  In  reality, 
it  was  found  that  comparing  the  difference  in  displacement 
for  the  system  effectiveness  resulted  in  the  error  value 
being  overly  sensitive  to  the  planfcrm  orientation  and 
relatively  insensitive  to  the  planform  curvature.  Each  sub¬ 
case  was  equally  weighted  in  the  summation  of  the  separate 
subcases.  The  subcase  error  value  was  a  function  of  the  dif 
ference  between  the  analytical  displacement  and  the  experi¬ 
mental  displacement.  Mathematically,  this  was  expressed 

as 

r  25 

NAS  ”  dij  EXP ^  2J  ^ 


z . 
J 


r 


Where  d(a,x)^j  NAg  represented  the  analytical  (NASTRAN)  dis¬ 
placements  which  were  functions  of  the  decision  variables  (a) 
and  the  state  variables  (x) ,  and  d^j  represented  the 

experimental  displacements.  The  analytical  displacements 
were  replaced  by  the  inverted  stiffness  matrix  and  forces 


d(a.£>ij  NAS 


I  S(a,x)ik  Pkj  „AS 


Where  the  stiffness  matrix  (S)  would  contain  the  decision 
variables  (a)  and  state  variables  (x)  .  The  load  cell 
spring  constant  was  included  as  a  decision  variable  in 
the  stiffness  matrix.  The  subcase  performance  indices 
can  be  reposed  mathematically  as 


r  25  r  25  -i  1 

j  =  Li=iL(k=l  Pkj  NAS)  '  dij  exp]2]2 


The  subcase  error  value  can  then  be  evaluated  for  a  solution 
set  by  varying  the  decision  variables,  subject  to  any  con¬ 
straining  equations  and  state  variable  relationships. 
However,  no  state  variable  relationships  were  identified  for 
the  problem  herein  and  no  constraints  were  imposed  in  the 
initial  searches  for  candidate  solutions.  Thus,  the  result¬ 
ing  error  function  was 


E  = 


10  rr  25  r  25  , 

Z  I  (  L  S(a  ,x)  •  r 

j  =  l  LLi=i  L  k=l  llc 


P  .  )  -  d .  . 

rkj  NAS''  ij  EXP 
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A  candidate  solution  set  was  used  in  the  reduced  error  value 
function  in  determining  the  system  effectiveness  for  the  deci¬ 
sion  variables.  The  error  calculated  represented  the  sum  of 
the  difference  in  calculated  vertical  displacements  versus 
experimental  displacements  for  250  data  points  (10  subcases, 

25  data  points  per  subcase).  The  Alter  routine  used  in  the 
NASTRAN  analysis  for  compiling  the  comparison  data  (error  value) 
is  presented  in  Appendix  E.  The  error  for  an  individual  sub¬ 
case  was  useful  in  distinguishing  relative  error  between  span- 
wise  bending  and  chordwise  bending  (psuedo-torsion) .  The 
boundary  parameters  were  a  significant  factor  in  reducing  the 
error  value. 


Variable  Parameters  Investigated.  The  problem  of  tuning 
the  FEM  to  predict  the  experimental  results  was  segregated 
into  two  distinct  problem  areas.  The  first  was  to  resolve  the 
dissimilarities  in  boundary  parameters  of  the  experimental  re¬ 
sults  and  the  analytical  results.  The  spring  constant  calcu¬ 
lated  for  the  actuator  arm  rotation  was  calculated  to  be  15.55 
x  10^  in  lb/rad.  This  was  the  result  of  the  actuator  arm 
length  of  6.0  inch  with  an  equivalent  spring  constant  of  3.2 
x  10^  lb/in  at  the  actuator  attachment  end  (Ref  9)  combined 
with  a  reaction  load  cell  spring  constant  of  5.0  x  10^  lb/in. 
However,  the  spring  constant  which  yielded  the  smallest  error 
value  (all  other  variables  held  constant)  was  1.15  x  10^  in  lb/ 
rad.  The  decision  to  use  the  boundary  parameter  values  which 
optimized  the  analytical  model  was  based  on  the  goal  of  op¬ 
timizing  the  planform  deformation  for  flutter  prediction  usage 
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and  the  experimental  boundary  conditions  were  not  those  of 
the  aircraft  configuration.  Once  an  optimum  set  of  boundary 
parameter  values  were  found,  then  the  deformation  of  the 
stabilator  planform  was  compared  to  the  experimental  dis¬ 
placements  for  optimization. 

The  variables  used  to  establish  the  boundary  conditions 
for  rigid  planform  tilt  and  rotation  were  bending  inertia  of 
the  torque  tube  sections  (part  of  the  hinge  fitting)  and  the 
spring  constant,  respectively.  They  were  iteratively  ad¬ 
justed  as  the  planform  variables  were  altered.  Three  basic 
parameters  were  investigated  with  respect  to  optimizing  the 
FEM  planform  for  analytical  displacement  prediction.  They 
were  inertias  of  the  membrane  elements  in  the  spanwise  and 
chordwise  directions  and  the  ratio  of  thickness  to  chord 
length.  Bending  in  the  chordwise  direction  was  associated 
with  a  psuedo-torsion  of  the  airfoil  to  distinguish  it  from 
bending  spanwise. 

Previous  investigations  were  accomplished  using  iso¬ 
tropic  elements.  Appendix  F  develops  the  theory  illustrating 
the  need  for  orthotropic  material  properties  to  model  the 
torsional  shear  flow  when  the  FEM  becomes  a  two-dimensional 
(planform)  model  (Fig  9).  Computations  were  made  using  iso¬ 
tropic  and  orthotropic  material  properties. 

The  thickness  to  chord  ratio  variable  was  used  as  an 
effective  weighting  factor  for  proportioning  both  bending 
and  psuedo-torsion  effects.  Since  the  thickness  was  di¬ 
vided  by  the  chord  length,  any  increase  to  the  thickness 


37 


with  Torsion  Cell 


without  Torsion  Cells 


Figure  9.  Torsion  Cell  Versus 

Two-Dimensional  Model 

was  equivalent  to  changing  the  airfoil  shape.  As  an 
example,  an  increase  in  this  variable  resulted  in  pro¬ 
portional  increases  in  bending  inertias  (both  directions) 
at  the  thickest  (Y  direction)  point  on  the  chord,  tapering 
to  a  minimum  (zero)  at  the  edges  for  a  particular  chord. 
Since  the  airfoil  also  tapers  from  the  root  rib  to  the  tip 
rib,  the  effective  weighting  in  thickness  was  also  realized 
spanwise . 


IV  Results 

Although  the  stabilator  was  restrained  in  a  relatively 
rigid  steel  and  aluminum  structure,  some  displacement  error 
was  evident  from  the  hysteresis  recorded  during  the  experi¬ 
mental  tests.  The  data  hysteresis  was  zero  in  most  cases 
and  repeatability  was  excellent.  The  maximum  deviation  for 
any  particular  load,  displacement  combination  was  found  to 
be  .004  inch.  The  data  was  averaged  to  arrive  at  the  final 
results  presented  in  Appendix  A.  The  tabular  results  were 
then  used  as  base  line  static  displacement  data  for  the 
Series  3  stabilator  and  all  subsequent  results  herein.  Static 
analysis  of  the  reaction  load  at  the  load  cell  was  found  to 
be  in  excellent  agreement  with  the  applied  load. 

Many  computations  were  made  using  the  scalar  error  value 
as  a  measure  of  the  accuracy  of  the  FEM  for  static  displace¬ 
ment.  Since  the  error  value  represented  the  total  difference 
(in  inches)  between  NASTRAN  predicted  displacements  and  ex¬ 
perimental  displacements  of  the  planform,  it  was  found  that 
variations  in  the  boundary  elements  (torque  tube  and  load 
cell)  overwhelmed  the  effects  of  varying  the  physical  para¬ 
meters  of  the  planform.  There  was  a  problem  of  having  too 
many  subcases  to  match  displacements,  and  not  enough  physical 
relationships  to  equate  some  of  the  subcases  (reducing  the 
independent  comparisons),  thus,  implying  the  scalar  error 
function  was  too  sensitive  to  planform  orientation  and  in¬ 
sensitive  to  the  problem  of  optimizing  the  planform  deformation. 
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The  goal  was  to  determine  the  variables  and  their  mag¬ 
nitudes  for  optimizing  the  two-dimensional  FEM  for  eventual 
modal  vibration  analysis.  The  boundary  element  variables 
were  reevaluated  after  each  iteration  of  the  planform  vari¬ 
ables  to  determine  if  further  improvement  in  the  scalar 
error  value  could  be  made.  The  boundary  parameter  values 
producing  the  best  scalar  error  value  were  significantly 
different  than  those  modelling  available,  reliable  data. 

As  stated  earlier,  the  rotational  spring  located  on  the 
torque  tube  to  restrain  the  rotation  was  computed  as  15.55  x 
10^  in  lb/rad  but  the  lowest  scalar  error  value  was  found 
when  a  value  of  1.18  x  10^  in  lb/rad  was  used.  A  comparison 
of  the  displacements  at  grid  point  (GP  151)  was  used  to 
optimize  the  bending  inertia  necessary  in  the  torque  tube. 

It  was  found  that  the  torque  tube  needed  to  be  more  flexible 
in  bending  by  a  factor  of  approximately  0.8  for  the  best 
planform  orientation.  This  value  was  based  on  all  other 
variables  being  held  constant,  and  the  optimum  spring  con- 
s  tant . 

With  the  above  changes  in  boundary  variables,  the  plan- 
form  of  the  stabilator  was  investigated.  Variables  used 
to  optimize  the  planform  curvature  were  orthotropic  bending 
moduli  material  properties  for  spanwise  and  chordwise  changes 
and  the  thickness  to  chord  ratio.  Curvature  of  the  planform 
chordwise  and  spanwise  could  be  altered  independently  using 
the  plate  bending  moduli.  Thickness  to  chord  ratio  changes 
were  used  when  comparison  of  the  data  indicated  the  center 
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(approximately  along  the  spar)  should  be  stiffer  (more 
heavily  weighted)  than  the  edges.  FEM  displacement  errors 
at  the  four  corners  of  the  planform  and  several  middle  grid 
points  were  used  to  determine  the  bending  and  psuedo-torsion 
changes  needed.  Effects  of  thickness  to  chord  ratio  changes 
were  investigated  and  found  to  adversely  affect  the  results 
obtained  using  orthotropic  properties. 

The  overall  results  indicated  it  was  necessary  to  change 
the  spanwise  bending  modulus  by  increasing  its  value  by  two 
percent.  The  chordwise  bending  modulus  (psuedo-torsion) 
was  increased  by  200  percent  (or  3  times  the  original  value). 
The  ratio  of  the  thickness  to  chord  length  ratio  was  not 
changed.  When  combined  with  the  optimum  torque  tube  bending 
inertia  and  actuator  arm  spring  constant,  the  scalar  error 
was  found  to  be  1.6  (i.e.  1.6  inches  of  displacement  error 
in  250  data  points).  .Prior*  to  any  optimization  attempts,  a 
scalar  error  of  10.4  was  found.  Optimization  of  the  boundary 
elements  accounted  for  lowering  the  scalar  error  to  1.8, 
prior  to  investigating  the  planform  variables.  Graphic 
presentation  of  the  NAjISTRAN  analytical  results  is  presented 
in  Appendix  G. 
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V  Conclusions  and  Recommendations 


The  thesis  investigated  a  two-dimensional  FEM  for  the 
prediction  of  static  displacement  for  a  T-38  horizontal 
stabilator  using  NASTRAN.  Experimental  data  was  obtained 
for  a  Series  3  stabilator.  The  data  provided  an  adequate 
base  of  information  for  comparison  purposes.  A  method  of 
comparison  was  developed  following  current  optimization 
principles.  The  scalar  error  value  gave  a  measure  of  the 
effectiveness  of  the  changes  imposed  on  the  combination  of 
boundary  parameters  and  planform  parameters.  It  was  found 
that  the  error  function  was  too  sensitive  to  planform  orien¬ 
tation,  since  a  lower  scalar  error  value  was  not  an  indication 
that  all  subcases  were  being  optimized.  This  feature  was 
used  to  determine  the  optimum  boundary  element  values  by 
finding  the  boundary  element  values  which  minimized  all 
subcases . 

Orthotropic  material  properties  were  essential  to  the 
ability  of  tuning  the  model  bending  properties  in  the  span- 
wise  and  chordwise  directions  independently.  This  was  not 
used  in  previous  theses  efforts.  The  optimization  process 
considered  relatively  few  variables  for  use  in  tuning  the 
model,  considering  the  large  number  of  elements  and  combina¬ 
tions  which  were  used  developing  the  FEM.  However,  many 
computations  were  required  to  narrow  the  results  to  an 
apparent  optimum  solution  set.  Based  on  the  original  model 
analyzed  and  the  tuning  technique  used  (scalar  error  value 


optimization) ,  the  FEM  using  orthotropic  plate  elements 
(appropriately  proportioned)  represents  a  good  stabilator 
finite  element  model. 

One  data  point  not  acquired  during  the  experimental 
tests  was  planform  rotation.  This  would  have  been  of 
great  benefit  in  determining  the  optimum  boundary  para¬ 
meter  value  for  establishing  planform  rotation.  It  is 
recommended  any  further  tests,  which  parallel  the  tests 
herein,  include  a  method  for  recording  the  planform  rota¬ 
tion  . 

It  is  recommended  the  results  of  this  thesis  be  con¬ 
sidered  for  vibration  analysis  of  the  T-38  stabilator. 

It  is  recommended  that  further  investigation  of  the  torque 
tube  (part  of  the  hinge  fitting)  be  performed  for  its 
contribution  to  stabilator  flutter.  This  recommendation 
is  based  on  the  large  relative  changes  in  the  scalar  error 
value  attributed  to  varying  the  torque  tube  bending  inertia 
and  the  attached  load  cell  spring  constant. 
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Appendix  A 

Experimental  Test  Results 

The  following  tables  present  the  vertical  deflections 
of  the  stabilator  due  to  static  loading.  The  location 
identification  is  given  in  terms  of  the  NASTRAN  grid  point. 
Reference  is  also  made  to  Northrop  static  displacement  test 
locations  in  parentheses,  which  were  used  to  establish  load 
and  deflection  points. 

The  maximum  applied  load  is  presented.  The  applied  load 
was  the  maximum  recorded  load  minus  the  weight  of  the  pad 
(8x8x1  inch  steel  plus  rubber  and  attaching  hardware, 
approximately  18  pounds).  The  reaction  load  was  that 
measured  by  the  load  cell  at  the  actuator  arm  for  the 
maximum  applied  load. 

The  deflection  data  represents  the  average  of  the 
vertical  displacements  found  for  successive  maximum  loads. 
Positive  is  upward  displacement,  negative  is  down. 
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TABLE  IV 


’001 


TABLE  IV  (cont) 


TABLE  IV  (cont) 


¥  V 


Appendix  B 

Graphic  Presentation  of  Experimental  Test  Results 

A  separate  FORTRAN  program  was  used  to  generate  a  dis¬ 
placement  deck  which  could  be  read  by  GCSNAST  (a  NASTRAN 
graphics  support  code) .  The  program  expanded  the  experi¬ 
mental  displacement  mesh  size  from  twenty-five  nodes  to 
166  grid  points  using  linear  interpolation  between  the  nodes. 
The  output  results  of  the  program  were  then  used  as  input  to 
GCSNAST  for  graphic  presentation  of  the  experimental  results 
presented  in  this  appendix. 

The  views  illustrated  present  the  same  view  as  the  ex¬ 
perimental  tests  set-up  (shown  in  Figs  2  and  3).  The  two- 
dimensional  FEM  forms  a  reference  plane  (horizontal  line) 
for  the  undeformed  stabilator,  such  that  the  leading  edge  is 
to  the  left  and  the  trailing  edge  is  to  the  right.  Reference 
loading  is  maximum  recorded  load  (not  applied  load  value). 
Load  position  is  also  noted. 
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EXPERIMENTAL  DEFORMATION  DUE  TO  A  4»l .  LI  LOAD  AT 
NODE  3  CHIC  POINT  14?  ul A  LINEAR  INTERPOLATION 


Figure  B-l.  Subcase  1  Experimental  Deformation 


EXPERIMENTAL  DEFORMATION  DUE  TO  A  4««.  LI  LOAD  AT 
NODE  ?  /  OfilD  POINT  IS1  UlA  LINEAR  INTERPOLATION 


igure  B-2.  Subcase  2  Experimental  Deformation 


EXPERIMENTAL  DEFORMATION  DUE  TO  A  78®  L)  LOAD  AT 
NODE  8  GRID  POINT  IS?  UIA  LINEAR  INTERPOLATION 


case 


EXPERIMENTAL  DEFORMATION  due  to  a  4«e  II  LOAD  AT 
NODE  9  /  GRID  POINT  163  OIA  LINEAR  INTERPOLATION 


Experimental  Deformation 


£KAt«  I  MENTAL  DEPORNATION  DUE  TO  AM  ll«  LI  LOAD  AT 
NODE  12  '  GRID  POINT  1SE  UIA  LINEAR  INTERPOLATION 


Subcase  5  Experimental  Deformation 


EXPERIMENTAL  DEFORMATION  DUE  TO  a  12M  II  LOAD  AT 
NODE  1)  '  GRID  POINT  157  VIA  LINEAR  INTERPOLATION 


EHPER1HCNTAL  OEPORNATION  DUE  TO  AN  844  LI  LOAD  AT 
MODE  14  -  GRID  POINT  1SI  UlA  LINEAR  INTERPOLATION 


Subcase  7  Experimental  Deformation 


EXPERIMENTAL  DEFORMATION  DUE  TO  A  1£«»  LI  LOAD  AT 
NODE  iy  '  GRID  POINT  160  UIA  LINEAR  INTERPOLATION 


Figure  B-8.  Subcase  8  Experimental  Deformat 


EXPERIMENTAL  DEFORMATION  DUE  TO  A  ISO*  L»  LOAD  AT 
MODE  18  '  GRID  POINT  184  UIA  LINEAR  INTERPOLATION 


Figure  B-9 .  Subcase  9  Experimental  Deformation 


EXPERIMENTAL  DEFORMATION  DUE  TO  A  I2H  t*  LOAD  AT 
NODE  19  /  OR  ID  POINT  t<6  UtA  LINEAR  INTERPOLATION 


Figure  B-IO.  Subcase  10  Experimental  Deformation 


As  mentioned  in  the  text,  a  separate  FORTRAN  program 
(BDATA)  was  used  to  generate  much  of  the  FEM.  The  program 
had  been  authored  by  the  thesis  advisor  (Capt  H.C.  Briggs). 

It  had  also  been  modified  extensively  during  its  usage  by 
subsequent  thesis  students,  Lassiter  and  Thomson.  The  pro¬ 
gram  generated  the  cards  defining  the  physical  properties  of 
the  stabilator  based  on  information  contained  in  Northrop 
Reports  (Ref  4  and  9)  and  Air  Force  Documents  (Ref  10)  for 
use  in  the  Bulk  Data  Deck.  Separate  algorithms  were  used  in 
forming  modules  to  compute  or  write  specific  segments  for 
output . 

The  modules  (or  sections)  used  in  the  program  in  order 
of  occurrence  can  be  described  as  follows: 

Module  1  -  Input  of  data  arrays  for  dimensions  and  function 

statements  for  other  algorithms. 

Module  2  -  User  decision  parameters  for  type  of  output 

(card-punch,  hard  copy-print,  or  tape),  type 
stabilator  mesh  (Series  3,  influence  coeffi¬ 
cient),  size  mesh  (number  of  elements  chordwise 
and  spanwise),  and  variables  used  for  refining 
displacement  (torsion,  thickness,  and  bending). 

Module  3  -  Grid  point  mesh  and  coordinates  are  calculated 


▼ 


for  following  algorithms.  Root  rib  coordinates 
are  found  from  interpolation  which  accounts  for 
the  fuselage  moldline.  The  influence  coefficient 
nodes  are  not  calculated,  but  specific  values 
established  in  data  tables. 


Module  4 


Module  5 


Grid  points  are  sequenced  for  CQUAD1  elements 
with  CTRIA1  elements  replacing  those  necessary 
to  construct  the  influence  coefficient  model. 
Connectivity  and  property  cards  are  developed 
for  CQUAD1  and  CTRIA1  elements. 

Grid  point  connectivity  for  CBAR  elements  for 
the  main  spar  and  associated  property  cards. 
These  are  explicit  for  the  influence  coeffi¬ 
cient  model. 


Module  6 


Module  7 


Module  8 


Module  9 


Generate  CBAR  elements  for  the  tip  rib,  con¬ 
nectivity,  and  associated  property  cards. 

Generate  CBAR  elements  for  the  leading  edge  extru¬ 
sion,  connectivity,  and  associated  property  cards. 

Generate  CBAR  elements  for  the  auxiliary  spar, 
connectivity,  and  associated  property  cards. 

Generate  elements  for  the  ribs  on  Series  2  only 
stabilator  (not  used  herein,  for  Series  3  analysis) 


Module  10 


Generate  CBAR  elements  for  the  root  rib,  con¬ 
nectivity,  and  associated  property  cards. 
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Module  11 


Connectivity  and  property  cards  for  the  torque 
tube  are  output  explicitly  based  on  model 
chosen  (influence  coefficient  versus  variable 
size  mesh) . 

Module  12  -  Generate  CBAR  elements  for  trailing  edge  clo¬ 

sure,  connectivity,  and  associated  property 
cards . 

Module  13  -  Output  coordinate  system  cards.  Two  rectangular 

systems  are  used,  the  global  system  and  a  sys¬ 
tem  embedded  in  the  planform  of  the  stabilator 
(C0RD2R) . 

Module  14  -  Output  GRID  cards.  Influence  coefficient  model 

grid  points  are  moved  or  generated  as  necessary, 
based  on  fixed  mesh  size. 

Module  15  -  Output  MAT1  and  MAT2  material  property  cards. 

These  are  explicit  except  for  variations  allowed 
in  determining  properties  for  the  orthotropic 
material . 

Module  16  -  Subroutine  for  interpolation  used  in  above 

algorithms . 

Module  17  -  Subroutine  for  development  of  property  card 

specifics  used  in  above  algorithms. 

Module  18  -  Subroutine  for  output  of  formatted  data  used 

in  above  algorithms. 
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Although  many  modules  were  used  in  the  generation  of  the 
data,  the  program  modular  format  was  of  a  great  benefit  during 
analysis  for  creating  Bulk  Data  Decks  with  varying  parameters. 
The  program  was  modified  for  this  analysis  to  include  the 
generation  of  those  elements,  connectivities,  and  grid  points 
relevant  to  the  generation  of  the  influence  coefficient  model. 
The  influence  model  was  developed  for  an  8  element  chordwise 
by  15  element  spanwise  mesh.  In  reviewing  the  contents  of  the 
program,  a  significant  error  in  calculating  the  inertia  of  the 
membrane  elements  was  discovered  and  corrected. 
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Appendix  D 

Grid  Point  Identifications 


Expanded  elemental  sections  from  Figure  7,  page  19,  are 
presented  in  this  appendix  for  better  clarification  of  the 
relative  locations  of  the  various  elements  used  in  developing 
the  FEM.  In  all  views,  the  tip  is  located  at  the  top  and 
forward  is  to  the  left. 
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SERIES  3  Itif  CO  MODEL 

GRID  POINTS  FOR  HSS  77.36  TO  8S.44 


Figure  17  1.  Grid  Point  Identification  GP  1  to  GP  27 


SERIES  1  INF  CO  MODEL 

GRID  POINTS  FOR  MSS  <6.8®  TO  68.88 


Figure  D-3.  Grid  Point  Identification  GP  55  to  GP  99 


Figure  D-4,  Grid  Point  Identification  GP  91  to  GP  126 


ERICS  3  IMF  CO  riODEl 

RIO  POINTS  FOR  HSS  33.64  TO  3S .  34 


Figure  D~5.  Grid  Point  Identification  GP  118  to  GP  166 


Appendix  E 


Alter  Routine  for  Error  Value  Calculation 


NASTRAN  Executive  Control  Deck  Section 

ID  JO  SAWDY 
APP  DISPLACEMENT 
SOL  1 
TIME  4 
ALTER  134 

PARTN  UGV , , ROWPT  /  NAS,,,  /  C,N,1  $ 

ADD  NAS, EXP  /  DIFF  /  C , Y , ALPHA = ( 1 . 0 , 0 . 0 )  /  C , Y , BETA= ( - 1 . 0 , 0 . 0 )  $ 
TRNSP  DIFF  /  DIFT  $ 

MPYAD  DIFT, DIFF,  /  DSQ  /  C,N,D  $ 

DIAGONAL  DSQ  /  ERROR  /  C , Y , OPT=COLUMN  /  V,Y,POWER=.5  $ 

MPYAD  SUMCOL, ERROR,  /  ERRSUM  /  C,N,0  $ 

MATPRN  NAS, EXP, DSQ, ERROR, ERRSUM  //  $ 

ENDALTER 

CEND 


The  above  Executive  Control  Deck  section  was  used  in 
formulating  the  scalar  error  value  in  NASTRAN.  After  the 
solution  of  displacements  for  all  subcases  is  computed,  the 
program  was  altered  (after  step  134)  to  extract  (PARTN) 
the  twenty-five  data  points  for  each  subcase  from  the 
solution  set  (UGV)  using  a  column  matrix  of  ones  and  zeros 
(ROWPT).  Using  the  ADD  operation,  the  experimental 
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displacements  (EXP)  are  subtracted  from  the  appropriate 
NASTRAN  computed  displacements  (NAS)  to  form  a  matrix  of 
difference  in  displacements  (DIFF).  The  transpose  (TRNSP) 
of  DIFF  is  formed  as  DIFT.  The  product  of  DIFT  and  DIFF 
matrices  then  forms  the  sura  of  the  squares  (DSQ) .  The 
square  root  of  the  diagonal  terms  was  formed  into  a  column 
vector  (ERROR)  representing  the  error  value  for  each  of  the 
subcases.  Using  the  operator  multiply  and  add  (MPYAD) ,  the 
subcase  error  values  (ERROR)  were  summed  using  a  row  vector 
(SUMCOL)  to  form  the  overall  scalar  error  value  (ERRSUM). 
This  value  represented  the  sum  of  the  differences  in  dis¬ 
placement  for  250  data  points  compared. 


Appendix  F 


Two-Dimensional  Model  Versus  Torsion-Cell 


Inasmuch  as  the  two-dimensional  model  is  made  of  solid 
elements,  it  is  deficient  the  mechanism  for  modelling  torsion¬ 
cell  shear  flow  contribution  to  the  stiffness  of  the  airfoil. 

A  psuedo-tors ional  mechanism  can  be  introduced  to  the  FEM  by 
considering  the  plates  as  orthotropic  elements  instead  of 
isotropic  elements  (as  was  utilized  in  previous  theses).  The 
application  of  orthotropic  elements  for  a  unit  span  is  in¬ 
vestigated  . 

Allow  the  cross-section  of  the  airfoil  to  be  approxi¬ 
mated  by  a  diamond  shape  of  height  h  and  chord  width  C  with 
a  constant  skin  thickness  of  t  as  shown. 

H 

y 

X 

(z  axis  out) 

Figure  F-l.  Diamond  Shape  Approximation  of  Airfoil 

The  two-dimensional  model  is  similarly  approximated  by 
a  thin  rectangular  shape  of  height  t  and  chord  width  C. 


Figure  F-2 .  Thin  Rectangular  Shape  Approximation 
for  FEM 
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It  is  desired  to  construct  the  El  and  GJ  products  of  the 
thin  rectangular  plate  such  that  they  are  equivalent  to 
those  of  the  diamond  shaped  airfoil.  The  EIa  product  of 
the  airfoil  surface  plates  about  the  x  axis  can  be  approxi¬ 
mated  by 


El 


where  ha(x)  " 


4E 


/ 


5/2  t  (ho(x))2dx 


a '  a 


-h 


x  + 


(Ref  11) 

(0  <  x  <  C/2) 


Substituting  the  equation  for  the  skin  surface  geometry,  the 
EIa  product  can  be  reduced  to 


Similarly,  the  El  product  for  the  thin  rectangular  shape 
about  the  x  axis  can  be  expressed  by 


El 


(Ref  12) 


For  a  constant  E,  an  equivalent  thickness  can  be  found  from 
equating  EIa  and  El  . 

t  =  3J  t  h2 
P  If  a  a  _ 

2 

Thus,  the  thickness  of  the  thin  rectangular  shape  is  a  func¬ 
tion  of  the  skin  thickness  and  airfoil  height.  This  relation¬ 
ship  is  used  in  developing  a  consistent  torsional  moment  of 
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inertia  for  G  J  and  G  J  products  about  the  z  axis.  The 

cL  ci  p  D 

airfoil  GaJa  product  can  be  expressed  by 


G  J 


a  a 


(s  is  airfoil  perimeter) 
(Ref  11) 


Substituting  for  the  airfoil  area  and  perimeter. 


4(c  -r-); 


Then  for  h  <<C,  the  product  G  J  reduces  to 

a  a  a 


The  GpJp  product  for  a  thin  rectangular  plate  can  be  ex¬ 
pressed  by 


GnJn 

P  P 


(Ref  12) 


Substituting  for  t  in 
equating  GflJa  to  GpJp, 


terms  of  the  airfoil  parameters  and 
an  equivalent  G  relationship  can  be 


deve loped . 


G  =  3Ga 
P  a 

Thus,  the  plate  shear  modulus  requires  three  times  the 
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the  airfoil  shear  modulus  value  for  an  equivalent  thickness 
FEM.  The  two-dimensional  model  did  not  use  G  (shear  modulus) 
for  the  plate  elements.  Therefore,  Young's  modulus  (E)  was 
varied  for  changes  in  bending.  However,  varying  E  for  an 
isotropic  plate  would  result  in  bending  changes  for  both 
directions  of  the  planform.  Thus,  it  became  imperative  to 
segregate  the  two  directions  of  bending  for  the  two-dimen¬ 
sional  FEM.  Orthotropic  plates  allowed  the  selection  of 
two  separate  E  values  for  optimizing  the  planform  deforma¬ 
tion. 
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Graphic  Presentation  of  NASTRAN  Analytical  Results 


The  views  presented  were  developed  using  the  output  of 
the  NASTRAN  analytical  displacements  for  all  (166)  grid  points 
as  input  to  the  graphics  code  GCSNAST.  The  orientation  of 
the  view  is  the  same  as  that  in  Appendix  B  to  facilitate  visual 
comparisons.  As  in  Appendix  B,  the  view  is  the  same  as  the 
experimental  tests  set-up  (shown  in  Figs  2  and  3).  The  two- 
dimensional  FEM  forms  a  reference  plane  (horizontal  line)  for 
the  undeformed  stabilator.  The  leading  edge  is  to  the  left 
and  the  trailing  edge  is  to  the  right. 


78 


Figure  G-l.  Subcase  1  NASTRAN  Analytical  Deformation  Due  to 
a  382  lb  Load  at  GP  147 


Figure  G-3.  Subcase  3  NASTRAN  Analytical  Deformation  Due  to 
a  682  lb  Load  at  CP  152 


igiire  G-4 .  Subcase  4  NASTRAN  Analytical  Oeiorination  Due 
a  382  lb  Load  at  CP  133 


Figure  C-5.  Subcase  5  NASTRAN  Analytical  Deformation  Due  Lo 
a  782  lb  Load  at  GP  156 
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